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Abstract. We calculate the quartet potential energy surface for Li+NH and use it to calculate elastic and 
spin-relaxation cross sections for collisions in magnetically trappable spin-stretched states. The potential 
is strongly anisotropic but spin-relaxation collisions are still suppressed by centrifugal barriers when both 
species are in spin-stretched states. In the ultracold regime, both the elastic and inelastic cross sections 
fluctuate dramatically as the potential is varied because of Feshbach resonances. The potential-dependence 
is considerably reduced at higher energies. The major effect of using an unconverged basis set in the scatter- 
ing calculations is to shift the resonances without changing their general behaviour. We have calculated the 
ratio of elastic and spin-relaxation cross sections, as a function of collision energy and magnetic field, for a 
variety of potential energy surfaces. Most of the surfaces produce ratios that are favorable for sympathetic 
cooling, at temperatures below about 20 mK. 

PACS. 34.50. Cx 



1 Introduction 

Since 1995, when the first gaseous Bose-Einstein conden- 
sates (BECs) were created Q3Q], the ability to cool, trap, 
and control atoms at ultracold (sub-millikelvin) tempera- 
tures has increased dramatically. More recently, attention 
has turned to cold and ultracold molecules J3j[4l[5] . With 
their more complex energy structure, ultracold molecules 
open up many new possibilities in fields such as high- 
precision measurement 0(7] , quantum computing [5J , and 
ultracold chemistry [5]. However, this extra complexity 
comes at a price, with the extra structure making molecules 
more difficult to cool. 

Atoms are usually cooled to microkelvin temperatures 
by laser Dopplcr cooling [10]. However laser cooling is 
much less effective for molecules than for atoms. For this 
reason, many other methods for cooling molecules have 
been developed. These can be categorized into two dis- 
tinct types, direct and indirect methods [TTlll2j . Indirect 
methods involve the formation of ultracold molecules from 
a sample of prccoolcd ultracold atoms, and deeply bound 
and absolute ground-state alkali-metal dimers have re- 
cently been formed, by magnetoassociation (Feshbach res- 
onance tuning) followed by laser state transfer with Stimu- 
lated Raman Adiabatic Passage (STIRAP) [T5HH1IMI5] ■ 
Indirect methods are limited by the range of ultracold 
atoms available, mainly to molecules formed from alkali- 
metal atoms. Direct methods, such as buffer-gas cooling 
[17] and Stark deceleration [11], involve directly cooling 
a sample of warm molecules down towards ultracold tem- 
peratures. NH molecules have been buffer-gas cooled and 



magnetically trapped in their X 3 I7~ state, both alone 
18,19,20 and in combination with N atoms [3T], and 
Stark-decelerated and electrostatically trapped in their 
metastable a 1 A state [35]. Proposals exist to transfer the 
metastable a 1 A molecules into the X 3 S" ground state, ei- 
ther with an intense laser field [33] or by collisions with Rb 
atoms [24] , However, the lowest temperatures currently 
obtainable by direct methods are tens to hundreds of mil- 
likelvin, so that a secondary cooling stage is required to 
reach ultracold temperatures. 

One possible secondary cooling method is sympathetic 
cooling, in which prccooled cold molecules are placed in 
contact with an ultracold atomic coolant and allowed to 
thcrmalizc. However, for sympathetic cooling to be effec- 
tive, the rate of elastic (cooling) collisions must greatly 
exceed the rate of inelastic (trap-loss and heating) colli- 
sions [25]. As a rule of thumb, the ratio of elastic to in- 
elastic cross sections must exceed 100. Molecules trapped 
in magnetic and electrostatic traps are held in low-ficld- 
seeking states (states in which the energy increases with 
increasing field) . Collisional spin- relaxation processes from 
low-field to untrapped high-field-seeking states will there- 
fore limit the efficiency of sympathetic cooling, and in 
many of the systems theoretically studied the rate of col- 
lisional trap loss will be too rapid for sympathetic cooling 
to work. Sympathetic cooling has yet to be demonstrated 
for neutral molecular systems. However it has been used 
to cool trapped ions [261I27] , to produce overlapping Bose- 
Einstein condensates [55], and to cool fermionic [2H] and 
bosonic [30] atoms to quantum degeneracy 
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The first to consider the sympathetic cooling of molecul 
with alkali- metal atoms were Soldan and Hutson |31j , who 
examined the interaction between rubidium atoms and 
NH(X 3 Z' _ ). They found that the dispersion-bound states 
of RbNH that correspond to the Rb( 2 S)+NH(X 3 i;-) dis- 
sociation threshold are crossed by deeper ion-pair states, 
which correspond to the Rb + ( 1 S') + NH~( 2 77) dissoci- 
ation threshold. These ion-pair states introduce mecha- 
nisms for inelastic collisions and three-body recombina- 
tion that may hamper sympathetic cooling. Lara et al. 
subsequently investigated the Rb( 2 S , )+OH( 2 7T 3 / 2 ) system 
[32,33] and also found ion-pair states that crossed the 
RbOH covalent states. They also carried out full quan- 
tum calculations of field-free low-energy collisions on the 
coupled potential energy surfaces and found fast spin re- 
laxation collisions that would prevent sympathetic cooling 
for most states. More recently, detailed potential energy 
surfaces and field-free scattering calculations for Rb+NH 
and Cs+NH were performed by Tacconi et al. }3~4"ll35ll3l)] . 
who found large rotational inelasticities. 

Sympathetic cooling of NH3 and ND3 has also been 
considered. Zuchowski and Hutson |37| surveyed the in- 
teraction potentials of NH3 with alkali- metal and alkaline- 
earth atoms. They found that all the systems studied 
exhibited large anisotropics that would produce strong 
rotation-tunnelling inelasticity, so that sympathetic cool- 
ing is unlikely to be successful for NH3 in low-field-seeking 
states. Zuchowski and Hutson [35] went on to examine the 
collisions of NH3 and ND3 with Rb, showing that there is a 
good prospect of sympathetic cooling of ND3 in high- field- 
seeking states, even in collision with magnetically trapped 
atoms. 

Lara et al. |33j identified several desirable features for 
sympathetic cooling of molecules that are not in their 
absolute ground state. First, it is desirable to use light 
atomic cooling partners, which produce high centrifugal 
barriers that can suppress inelastic loss in channels with 
a small kinetic energy release. Secondly, it is desirable to 
use molecules where the electron spin is weakly coupled 
to the intermolecular axis, such as molecules with Hund's 
case (b) coupling. Thirdly, it is desirable to use atom- 
molecule systems for which the interaction potentials are 
only weakly anisotropic. They suggested that closed-shell 
atoms such as the alkaline earths might make good colli- 
sion partners for magnetically trapped molecules. 

In 2009 Soldan, Zuchowski, and Hutson [39] followed 
up these suggestions: they chose the case (b) molecule 
NH(X 3 Z"~) as a test case and investigated its interactions 
with both alkali-metal and alkaline-earth atoms. All the 
alkali-metal and the heavier alkaline-earth systems (Ca- 
NH and Sr-NH) were found to have strongly anisotropic 
interaction potentials, but Be-NH and Mg-NH were found 
to have weak anisotropy. All the systems have ion-pair 
states, but for Be and Mg the ion-pair states cross the 
dispersion-bound state high up on the repulsive wall, so 
that the ion-pair states will be energetically inaccessi- 
ble in low-temperature collisions. Wallis and Hutson [40] 
then performed detailed quantum scattering calculations 
on Mg+NH(X 3 I7 _ ) and showed that the ratio of elastic 



to inelastic cross sections is large enough to allow sympa- 
thetic cooling over a wide range of collision energies and 
magnetic fields. Sympathetic cooling of NH with Mg has a 
good prospect of success if the molecules can be precooled 
to 10 mK and good trap overlap can be achieved. 

Part of the success of Mg as a coolant comes from 
the fact that its ground state has zero electron spin, so 
that only spin relaxation and not spin exchange is pos- 
sible in collisions with NH. However, the lack of electron 
spin also prevents magnetic trapping. In addition, Mg is 
quite hard to laser-cool, though Mchlstaublcr et al. [4"T] 
have been successful in achieving sub-Doppler tempera- 
tures. An alternative way to prevent spin exchange is to 
use an alkali-metal atom as the collision partner, while en- 
suring that all collisions are between atoms and molecules 
in spin-stretched states, with the maximum possible val- 
ues of both electron and nuclear spin projection quantum 
numbers. The use of spin-stretched states also prevents 
strong interactions involving ion-pair states. It is impor- 
tant to know whether this approach would result in low 
enough inelastic collision rates, despite the much larger 
anisotropics in alkali-NH systems. Previous collision cal- 
culations on Rb+NH and Cs+NH [SB] did not include the 
effects of magnetic fields and did not consider the spin- 
changing collisions that are important for sympathetic 
cooling. 

In this paper we assess the prospects for sympathetic 
cooling of cold NH(X 3 Z ,_ ) molecules with ultracold Li( 2 S') 
atoms. We choose Li as the collision partner because its 
low mass produces high centrifugal barriers that suppress 
spin-changing collisions at low energies. 7 Li has nuclear 
spin /Li = 3/2, and has a magnetically trappable spin- 
stretched state with total angular momentum F = 2 and 
Mp = 2. The 14 NH molecule has complicated hyperfinc 
structure due to its electron spin 5 = 1 and two nu- 
clear spins In = 1 and In = 1/2; its magnetically trap- 
pable spin-stretched state has F = 5/2 and Mp = +5/2. 
Collisions between spin-stretched Li and NH have Ms = 
m su + m SNH — +3/2 and thus occur entirely on the quar- 
tet surface of Li-NH. We confine our attention to these 
collisions. Collisions between states with \Ms\ ^ 3/2 in- 
volve the doublet surface as well as the quartet surface 
and are likely to lead to fast spin exchange. 

The structure of this paper is as follows. Section [5] 
describes electronic structure calculations of the quartet 
potential energy surface for Li( 2 S')+NH(X 3 Z ,_ ). Section[3] 
describes quantum scattering calculations and Section 2] 
explores the prospects for sympathetic cooling. 

2 Potential Energy Surface 

We have calculated the quartet interaction potential for 
Li( 2 S')+NH(X 3 Z' _ ), using high-level ab initio electronic 
structure calculations. Calculations were performed using 
an open-shell spin-restricted version [52] of the coupled 
cluster method |43j with single, double, and perturbative 
triple excitations [RCCSD(T)]. We used the aug-cc-pVTZ 
basis set for N and H, and the uncontracted cc-pVTZ ba- 
sis for the Li atom [44) , augmented by Gaussian functions 
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Fig. 1. Li( 2 5)+NH(X 3 X'~) quartet potential energy surface. 
Contours are labeled in cm -1 . 6 = 0° corresponds to the linear 
N-H-Li geometry. 



with exponents s:0.001, p:0.008, d:0.034, and f:0.073. To 
improve the basis set for the Van der Waals dispersion 
interaction we included spd bond functions (with expo- 
nents 0.9, 0.3, 0.1 for s and p and 0.6 and 0.2 for d and 
f) at the Li-NH center of mass. All ab initio calculations 
were performed using the MOLPRO package [55] and cor- 
rected for basis-set superposition error using the counter- 
poise method of Boys and Bernardi . 

Calculations were performed on a grid in Jacobi coor- 
dinates (R, 9), with a fixed NH bond length r = 1.0367 A 
[47] . R is the separation between the NH center of mass 
and the Li atom, and 9 is the angle between the NH bond 
vector and R, with 9 = corresponding to the linear NH- 
Li geometry. The angular grid is a set of 11 angles cor- 
responding to Gauss-Lobatto quadrature points, and the 
radial grid extends from high up on the repulsive potential 
wall to 15 A. 

For each angle, radial interpolation for R < 15 A and 
extrapolation for R > 15 A is performed using the repro- 
ducing kernel Hilbert space (RKHS) method [4"51l4"9] to 
obtain the potential for any given value of R. The poten- 
tial is then expanded in Legendre polynomials, 



V(R,9) = J2Vx{R)Px{cos9), 



(1) 



for A < 8. The radial strength coefficients V\(R) are eval- 
uated by integrating over the angular grid at each R using 
Gauss-Lobatto quadrature. 

Figure [1] shows the calculated quartet potential energy 
surface for Li-NH. The global potential minimum is about 
1800 cm -1 deep and occurs at (9=180° (Li-NH). A sec- 
ondary minimum about 113 enr 1 deep occurs at 9=0°. 
The large anisotropy is due to s-p mixing of the Li orbitals, 



which is much stronger at the Li-NH geometry than for 
NH-Li [39]. 

Before examining the scattering dynamics, it is in- 
structive to consider some analytical properties of the po- 
tential. For a Van der Waals potential that decays at long 
range as —Cq/R 6 , a characteristic potential length may 
be defined l50l. 



p 1 f^C 6 

RvdW - 2 {— 



1/4 



(2) 



Alternatively, the average scattering length a, introduced 
by Gribakin and Flambaum [51) . is 



4tt 



r(i/4)' 



:R 



vdW 



0.956i?- 



vdW- 



(3) 



Gribakin and Flambaum showed that the scattering length 
follows 

1 - tan (<£--] , (4) 



S 



where 



<2> = 



[-2^V(R)/h 2 ] 1/2 dR 



(5) 



is evaluated at threshold, with Rq the inner classical turn- 
ing point. When a is large and positive, the energy of the 
highest bound state is just below threshold and is given 

by 

h 

2fi(a — ay 

If the potential is varied, by scaling or otherwise, the scat- 
tering length passes through a pole every time an s-wave 
bound state occurs at threshold. For Li( 2 S , )+NH(X 3 Z' _ ), 
C 6 = 164.9 E h a%, and thus i? vdw = 10.89 A and a =10.42 A, 
corresponding to an average elastic cross section of <7 c ias = 
Ana 2 = 1365 A 2 . 



3 Scattering Calculations 

The Hamiltonian for a 2 S atom colliding with a rigid-rotor 
3 S diatomic molecule may be written 



H 



h 2 d 2 



v- 



2/i dR 2 2[iR 2 



Hz + V SS (R) + V(R, 



(7) 

where L 2 is the space-fixed end-over-end rotation oper- 
ator for the atom and the diatomic molecule about one 
another, H m0 n is the field-free Hamiltonian for the di- 
atomic monomer, Hz represents the Zeeman interaction, 
V SS (R) is the anisotropic intcrmolccular spin-spin interac- 
tion, V(R, 9) is the intermolccular potential, and /i is the 
reduced mass of the colliding system. The approximation 
of treating NH as a rigid rotor is justified in more detail 
below. 

The field-free monomer Hamiltonian for NH(X 3 i7~) is 
H mon = b mt h 2 + 7, ls n • sn H 



+ A ss 



47T 

T 



\/6]r(-l) 9 y 2 - g (r) [s NH (8 SNH]?te) 
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\sa — sb I and sa + Sb- For Li-NH the allowed values are 
S = j (doublet) and § (quartet). The interaction po- 
tential Vi n t(-R, 9) may be written in terms of projection 
operators, 

V int (R,8)= J2 \S)V S (R,6){S\ (12) 



and the matrix element of V[ n t(R, 9) is 
(sAm sA SBm S Bmrn nB LM L \V in t(R,0)\sAm'sBm' sB n' B m' nB L'M' L } 

'nj - '2 

_ ^ V _^^2sa+2sb— TO e A- m,B- Ml _|_ 1) 



Fig. 2. 

function 

+ l,m SLi 



(color online). Li( 2 S')+NH(X 3 r~) thresholds as a 
of magnetic field. The spin-stretched (?n SNH = 
= +f) Ms = +f state is shown in red. 
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sa s B 



S 



m sA m ! 
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where n and snh are the NH rotation and electron spin 
angular momentum operators, 6 ro t = h 2 /2fj,r 2 non = 16.343 
cur 1 is the NH rotational constant, 7„ s = —0.055 cm" 1 is 
the spin-rotation constant and A ss = 0.92 cm -1 is the spin- 
spin interaction constant. The values of the spectroscopic 
constants for the vibrational ground state are from rcf. 
[521. The Zeeman Hamiltonian is 



Hz = 9e^B ■ (s N h + S L i), 



(9) 



where B is the magnetic field vector, fi B is the Bohr mag- 
neton and g e is the electron g-factor. The anisotropic in- 
termolecular spin-spin interaction can be expressed as 

V SS (R) = X(R) [s L i ■ s N H - 3(s Li • e R )(s N H ■ e R )] , (10) 

where Su and Snh are the spin angular momentum oper- 
ators for the Li atom and the NH molecule respectively, 
e R is a unit vector along R, and \(R) = —Eha 2 /(R/cto) 3 , 
where a ~ 1/137 is the fine-structure constant. 

We construct the collision Hamiltonian in the fully un- 
coupled basis set \nm n }\sjq B m s - NH }\su'msi ji )\LMi J ) , where 
rrii is the projection of the angular momentum term s^ 
onto the space-fixed magnetic field axis. The matrix el- 
ements for NH monomer operators are the same as for 
scattering of NH from a closed-shell atom [53], with the 
addition of factors S, t 
interaction has matrix elements 



The doublet and quartet interaction potentials Vs(R,0) 
have the same long-range coefficients, so become degen- 
erate when the N atom and NH molecule are far enough 
apart that their valence shells do not overlap. In the present 
work we approximate the operator Vi n t (i?, 9) by taking 
Vs — V3/2 f° r both spin states. This approximation is 
reasonable because we wish to focus on collisions between 
atoms and molecules in spin-stretched states, and V1/2 has 
no matrix elements (diagonal or off-diagonal) involving 
such states. When this approximation is made, orthog- 
onality relations for the 3j symbols reduce Eq. [LJ] to a 
form diagonal both in m s A and m sB . The explicit expres- 
sion for (n B m nB LML\Vs(R, 9)\n B m' nB L' M' L ) is the same 
as for scattering of NH from a closed-shell atom [53] , with 
the addition of factors S msAm ' . 

In a magnetic field, the lowest Li-NH threshold (N = 
0, snh = IjSLi = h) splits into six Zeeman sub-levels, as 



shown in Figure [2] There is one Ms 



state, two 



degenerate Ms 



states, two degenerate Ms = 



m 3 Am' a • The intermolecular spin-spin 



states, and a single spin-stretched Ms = | state. 

In a magnetic field the projection of the total angular 
momentum M to t = m n + Ms + Ml is conserved. Thus in 
the rotational ground state N = 0, any change in Ms 
must be accompanied by a change in M^. For s-wavc 
scattering (L = 0) from the spin-stretched Ms = +§ 
state, this requires that L — — >• L > 2. Thus the 
dominant outgoing channels are I M s =±U = 2) and 



I Ms = — 5 , L = 4) , which have centrifugal barriers that 
(s A m sA s B m sB n B m nB LM L \V ss \s A m sA s B rn sB n B m nB L M L ^ U pp re ss inelastic scattering at low energies and magnetic 
= V3Q\(R)8 ' S ' (— i) SA + SB-m s A-msB_Mj fi e kis. The heights of the centrifugal barriers are approxi- 
™ B " B m " Bm " B .mately = (h 2 L(L + l)/n)§ (54C 6 ) - £, which for LiNH 

[s A (s A + l)(2s A + l)s B (s B + l)(2s B + 1)(2L + l)(2L' + l)] 2 are E 2 f = 60 mK and £ c 4 f = 368 mK. 

2 L' \ For magnetically trapped spin-stretched states, colli- 
-Qi — 92 M' Jsional spin relaxation to other Zeeman states will result 
in heating and/or loss of molecules from the trap. Spin 
s B 1 Sb \ relaxation occurs via two mechanisms. The first is via the 
-m S B 92 m'sB 7 ai^cct coupling of the incident state to outgoing states 
through the anisotropic spin-spin interaction, by a similar 
The total spin S of a system made up of an open-shell mechanism to spin relaxation in collisions between alkali- 
atom and an open-shell molecule can take values between metal atoms. The second occurs via the interplay of the 
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these avoided crossings are at a distance -Rax such that 
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Fig. 3. (color online). S-wave total inelastic cross sections cal- 
culated as a function of collision energy for Rax > Rmax = 50 
A (blue, dashed) and Rax < Rmax = 600 A (red, solid) for 
magnetic fields of 1 and 10 Gauss. Rax(1 G) = 476 A and 
Rax(10 G) = 150 A. 



potential anisotropy and the intramolecular NH spin-spin 
interaction |54j . which results in an inelastic cross section 
proportional to l/b 2 . ot [551155] . 

We have carried out low-energy quantum scattering 
calculations using the MOLSCAT package [57], as mod- 
ified to handle collisions in magnetic fields [53]. Sets of 
coupled-channel equations are constructed for each M tot 
in the uncoupled basis set described above. Solutions to 
the coupled-channel equations arc then propagated over 
a grid in R, using the hybrid log-derivative method of 
Alexander and Manolopoulos [55], which in the short- 
range region (i? m in < R < -Rmid) uses a fixed-step-size log- 
derivative propagator and in the long-range region (i? m id < 
R < -Rmax) a variable step-size Airy propagator. At -R max , 
the log-derivative solutions arc transformed into the basis 
set of asymptotic eigenfunctions and matched to scatter- 
ing boundary conditions to obtain the 5-matrix S Mtot . 
The scattering cross sections between levels i and / are 
then obtained from 



\SifS 



iJ°L t L f 



a iLifL f 



M tot LiL f 



(14) 



where ki is the wave vector in the incident channel, k 2 = 

The anisotropic intcrmolecular spin-spin interaction 
directly couples the incident s-wave channel Mg = +|, L = 
to the two Ms = +h,L = 2 channels and to the (—1, +5), 
— h , L — 2 channel. Since the L = 2 channels have bar- 
riers at long range, there are narrowly avoided crossings 
between the channel adiabats [59j at very long range when 
the kinetic energy release is small (at low magnetic field) . 
If the long-range behaviour in both the L = and L = 2 
channels is of the form 



g e (i B BAM s 



h 2 [LfjLf + V-LijLi + l) 
2fiR 2 AX 



(16) 



where g^ji^BAMs is the energy separation of the two 
adiabats due to the magnetic field B, Li and Lf are the 
values of L in the incoming and outgoing channels, and 



M. 



Ms^ To account correctly for the intcr- 



molecular spin-spin coupling, the scattering solutions need 
to be propagated out beyond the largest value of -Rax- 
Figure [3J shows the effect of placing R max < Rax and 
-Rax < -Rmax on the total inelastic s-wave cross sections 
as a function of collision energy. We consider only mag- 
netic fields relevant to sympathetic cooling down to 1 G, 
at which point the largest -Rax = 476 A, and thus the 
propagation parameters i? m m = 1.8 A, R mid = 12.5 A 
and -Rmax = 600 A were used throughout. 

Low-energy scattering depends strongly on the details 
of the potential energy surface. Thus when converging the 
basis set for the scattering calculations wc need consider 
the uncertainty in the interaction potential as well as con- 
vergence with respect to the variables n max and L max that 
govern the basis set size [60] . To explore the convergence 
with respect to the potential we introduce a scaling fac- 
tor A sc i that allows the interaction potential to be varied 
within its error bounds, 



V 



scaled 



(R,e) = x acl v(R,e). 



(17) 



V(R) 



h 2 L(L + l) 
2/xi? 2 



R 6 



0(R~ 7 ), 



(15) 



Since the main features of the cross sections are deter- 
mined by the scattering length and the positions of bound 
states that cause resonances, both of which are tuned by 
A sc i, scaling the potential in this way is sufficient to ex- 
plore the range of possible collision behaviour. This also 
justifies fixing the NH bond length r in our calculations: 
the collisions arc approximately vibrationally adiabatic, 
and allowing r to relax would simply deepen the poten- 
tial slightly and correspond approximately to a value of 
Asci > 1. 

Figure 2] shows the potential-dependence of the s-wave 
spin-stretched clastic scattering cross section for the three 
basis sets n max = 6,i max = 6 (645 channels), n max = 
6, Lmax = 8 (937 channels), and n max = 8, i max = 8 (1403 
channels) at a collision energy of 10 -6 K and a magnetic 
field of 10 G. Calculations with n max = 8, L max = 10 
(1887 channels) gave results essentially identical to those 
with n max = 8,-L max = 8. The computer time scales as 
the number of channels cubed. It may be seen that, for all 
three basis sets, the elastic cross sections fluctuate about 
^|he _avcrage cross section a e i as and the total inelastic cross 
sections have peaks at corresponding positions. The fluc- 
tuations (resonances) have quite similar form for different 
basis sets, but the positions of the resonances are shifted. 
This is similar to the result obtained by Janssen et al. [60] 
for NH-NH. 

It is clear that basis sets smaller than n max = 8, L max = 
8 do not give "converged" results for a specific interaction 
potential. However, the calculated interaction potential is 
itself probably accurate only to about ±5%. Even a fully 
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Fig. 4. (color online). Spin-stretched s-wave elastic (upper 
panel) and inelastic (lower panel) cross sections calculated 
for basis sets with n max = 6, Z/ ma x = 6 (green, dotted), 
n m ax = 6, Lmax = 8 (blue, dashed), and n max = 8, L max = 8 
(red, solid), as a function of the potential scaling factor A sc i at 
a collision energy of 10 -6 K and a magnetic field strength of 
10 Gauss. The average cross section <7 e ias = 4ira 2 — 1365 A 2 is 
also shown (black, heavy-dots). 



converged basis set would therefore not give reliable pre- 
dictions of numerical values for elastic (or inelastic) cross 
sections in the ultracold regime. Nevertheless, even an un- 
converged basis set reliably predicts the probability that 
the elastic and inelastic cross sections will lie in a partic- 
ular range of values, and thus the likelihood that sympa- 
thetic cooling can succeed. Until the actual cross sections 
are measured, we have to accept this uncertainty in the 
results, and until then there is no point in using very ex- 
pensive basis sets that aim for complete convergence. 

The converged n max = 8,L max = 8 basis is compu- 
tationally too expensive for calculations involving p-wave 
and higher partial-wave contributions, which are required 
at temperatures above 10 /kK. The smaller and more com- 
putationally manageable n max = 6, L max — 8 has there- 
fore been used in the remainder of this work. Figure [S] 
shows the dependence of the elastic and total inelastic 
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Fig. 5. (color online). Dependence of the elastic (red) and total 
inelastic (blue) cross sections on the potential scaling factor 
Asci at collision energies of 1 fiK (top), 10 mK (center) and 
100 mK (bottom) at a magnetic field strength of 10 G. The 
s-wave contribution to each cross section is shown as a dashed 
line (and is the only significant contribution at 10 -6 K). The 
average s-wave cross section Oeias = 47ra 2 = 1365 A 2 is shown 
as a dotted black line. The calculations use the basis set with 

^max — 6, Lm&x — 8. 
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cross sections with this basis set for a wider range of A sc i 
for collision energies of 1 /iK, 10 mK and 100 mK: it may 
be seen that the higher partial waves that contribute at 
higher energies "fill in" the deep dips that occur in the s- 
wave cross sections, but that additional structure appears 
at higher energy because of resonances in higher partial 
waves. 



4 Prospects for Sympathetic Cooling 

In a simple hard-sphere model of sympathetic cooling, ne- 
glecting inelastic collisions, the temperature relaxes to- 
wards equilibrium and reaches a 1/e point after (mi + 
7712) 2 j1m,\m,i collisions [25], where m\ and 11J2 are the 
masses of the two species. For sympathetic cooling to be 
successful we need the ratio of elastic to inelastic cross 
sections to be much larger than this. The general rule of 
thumb is that the ratio of the elastic to the total inelas- 
tic cross sections, 7 = Ccias/tfmcias, must be greater than 
about 100. 

Figure [5] shows the spin-stretched elastic and total in- 
elastic cross sections, incorporating all s, p, and d partial 
wave contributions (L = 0,1,2) for the original ab initio 
potentials with A sc i = 1 and for representative scaled po- 
tentials with Asci = 0.96 and 0.915, chosen to give a = 0.3a 
and a = a, respectively. The potential with A sc i = 1 has 
a = —12.5a, and thus an atypically large value of the elas- 
tic cross section. The spd basis set gives convergence of the 
partial- wave sum in cross sections for collision energies up 
to about 60 mK. At low magnetic fields the kinetic energy 
release (i.e. the Zeeman splitting between the incoming 
and outgoings channels) is less than the height of the bar- 
riers in the outgoing channels, so that the inelastic cross 
section is suppressed. As the magnetic field is decreased, 
the Zeeman splitting decreases, the centrifugal suppres- 
sion increases, and the inelastic cross section decreases. 
The dominant elastic processes are AL = 0, so the elastic 
cross section does not display significant field-dependence. 
It may also be seen from Figure [SJ that the ratio 7 is in 
excess of 100 for low magnetic fields and a wide range of 
collision energies up to about 30 mK. However, due to the 
presence of the intermolecular spin-spin interaction |59j . 
the reduction of the inelastic cross section at low fields is 
not as dramatic as in Mg+NH [3D]. 

To assess the prospect of sympathetic cooling over the 
range of possible potentials, Figure [7] shows contour plots 
of 7 as a function of the collision energy and magnetic field 
strength for a variety of A values, chosen to give a range of 
scattering lengths from —12.5a to +12. 5a. It may be seen 
that the majority of possible potentials give favourable 
values of 7 over a wide range of fields and energies. The 
unfavourable potential at the upper right is also atypical: 
it has a = 0.1a, so that the s-wave elastic scattering cross 
section is a factor of 100 below the typical value of Ana 2 . 
Such unfavourable behaviour might occur for any system, 
but is not very likely. 

In an unbiased magnetic quadrupole trap, with a field 
zero at the center, the molecules in the m SNH = 1 state at 
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Fig. 6. (color online). Li+NH spin-stretched elastic (solid, 
red) and total inelastic (solid, blue) cross sections as a func- 
tion of collision energy for various magnetic fields, for the 
potentials with A sc i = 1, 0.96 and 0.915, corresponding to 
a — —12.5a, 0.3a, a, respectively. The cross sections include all 
s, p, and d-wave (L = 0, 1, 2) contributions. The s-wave elastic 
(dashed, red) and total inelastic (dashed, blue) cross sections 
are also shown. 
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Fig. 7. (color online). Contour plots of the ratio 7 of the elastic to total inelastic cross sections including all s, p and d-wave 
contributions. Contour plots are shown for six different values of the potential scaling factor (clockwise from the top left): (a) 
A sc i = 1, a = -12.5a. (b) A sc i = 0.9915, a = -a, (c) X scl = 0.97, a = 0.1a, (d) A sc i = 0.96, a = 0.3a, (e) X scl = 0.915, a = a, and 
(f) A sc i — 1.00078, a = 12.5a, The solid red lines show the maximum field sampled by the spin-stretched state of the molecules, 
B — 6fcBTNH/peA'B and the dashed red line shows a similar boundary for Li atoms at TLi = 1 mK. 



a temperature Tnh will have a Boltzmann density distri- 
bution 

Po \ kb ?nh / 

At a given temperature on the energy axis of the contour 
plots in Figure [7] less than 1% of molecules will experience 
a magnetic field greater than B = 6fcB?NH/3eMB, which is 
shown as a solid red line. If the atomic cloud occupied the 
same volume as the molecular cloud, the large majority of 
collisions would occur at energy/field combinations below 
this line. However, at the start of sympathetic cooling the 
Li is likely to be substantially colder than the NH, prob- 
ably 1 mK or less. The Doppler limit for Li is 140 /j,K. 
The magnetic moment of 7 Li in its F = 2, Mp = 2 state 
is only half that of NH, but nevertheless it is likely that 
the Li cloud will be significantly smaller. Li-NH collisions 
will take place only where both species are present, near 
the centre of the trap where fields are low. This will re- 
duce the frequency of collisions but may nevertheless be 
advantageous because it will mean that collisions mostly 
take place in regions well below the solid red line. Over 
99% of Li atoms at a temperature Tli will occupy regions 
of the trap with fields less than B = 12/03X^1/ g e pB, and 
this boundary is shown in Figure [7] as a dashed red line 
for Li at 1 mK. For many potentials the inelasticity in 
this region remains low to temperatures above 100 mK at 
fields below 100 G. 

It can be seen from Figure [7] that the contour plots 
depend mostly on the magnitude of the scattering length 
and not on its sign, because the s-wave elastic cross sec- 



tion is proportional to the square of the scattering length. 
There are small deviations from this: for a w a there is 
a p-wave resonance around 50 mK, and for a ~ —a there 
is a d-wave resonance around 50 mK and 500 G. If the 
real potential is such that a is accidentally close to zero, 
such as for A sc i = 0.97, sympathetic cooling will not suc- 
ceed because of the small elastic cross section. However 
as \a\/a increases to "typical" values close to 1 (such as 
for Asci = 0.96), sympathetic cooling has a good chance 
of success. For all potentials with \a\/a > 1 it can be seen 
that, once sympathetic cooling has started and the sam- 
ple of molecules is cooled to sub-millikelvin temperatures, 
7 increases and the molecules become increasingly stable 
to spin-relaxation. The temperature at which sympathetic 
cooling begins to be effective depends quite strongly on the 
potential, but for some potentials it is as high as 100 mK. 

Changing the reduced mass of the collision system has 
an effect quite similar to scaling the potential [61] . Thus if 
the scattering length for one isotopic combination proves 
to be too small to allow sympathetic cooling to succeed, 
another combination may well be more successful. 

A further important consideration in a sympathetic 
cooling experiment concerns the absolute collision rate. As 
seen above, the elastic cross section in the ultracold regime 
fluctuates dramatically near Feshbach resonances but its 
"background" value is 47ra 2 , which for Li+NH is 1365 A 2 . 
The elastic cross section drops off at higher energy for 
some potentials, but in most cases it does not change by 
more than a factor of 3 from 47ra 2 (in either direction) 
at collision energies around 100 mK. At this temperature, 
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with a reasonable Li atom density of 10 10 cm -3 and an 
elastic cross section of 10 3 A, each molecule will experience 
a collision with a Li atom once every 0.8 ms. However, this 
time will be extended if only part of the molecular cloud 
overlaps with the atomic trap. In a magnetic quadrupolc 
trap with a molecular temperature much higher than the 
atomic temperature, the overlap between the two clouds 
will be quite poor. This might be alleviated by using a trap 
of Ioffc-Pritchard type, with a finite magnetic field at the 
bottom of the trap. This would have the beneficial side- 
effect of preventing Majorana losses [52111)3"] at the center 
of the trap. As cooling progresses, the molecular cloud will 
collapse towards the center of the trap, so the overlap with 
the atomic trap will improve. As this occurs, the bias field 
of the Ioffe-Pritchard trap can be reduced to increase the 
density and the corresponding collision rate. 



5 Conclusions 

We have calculated the quartet potential energy surface 
for NH interacting with Li and used it to calculate elas- 
tic and inelastic cross sections for collisions between cold 
NH and Li in magnetically trappable spin-stretched states. 
The potential energy surface is highly anisotropic, with 
a well depth about 1800 cm" 1 at the Li-NH geometry 
but only 113 cm -1 and the NH-Li geometry. The strong 
anisotropy drives spin-relaxation collisions that are con- 
siderably faster than in systems such as Mg-NH and He- 
NH, but are nevertheless suppressed by centrifugal barri- 
ers when both species are in spin-stretched states. 

We have explored the dependence of the collision prop- 
erties on the potential energy surface, or equivalently on 
the reduced mass of the colliding system. In the ultracold 
regime, both the elastic and inelastic (spin-relaxation) cross 
sections show dramatic fluctuations as the potential is var- 
ied, due to Feshbach resonances that occur as a bound 
state passes through zero energy. The potential-dependence 
is considerably reduced at higher energies, when higher 
partial waves contribute to the scattering. The major ef- 
fect of using an unconverged basis set in the scattering cal- 
culations is to shift the positions of the resonances without 
changing their general behaviour. 

We have calculated the ratio of elastic and inelastic 
(spin-relaxation) cross sections, as a function of collision 
energy and magnetic field, for a variety of potential en- 
ergy surfaces. This ratio is key to the success of sympa- 
thetic cooling, because spin-relaxation cross sections con- 
vert internal energy into relative translation of the collid- 
ing species and thus cause trap loss. We find that most 
(but not all) potential energy surfaces produce ratios that 
are favorable for sympathetic cooling, provided the molecules 
can be precooled to a temperature around 20 mK. The 
exceptions are potentials where the scattering length is 
accidentally near zero, which produce small elastic cross 
sections in the ultracold regime. If this unlucky situation 
occurs in practice, it can probably be avoided by using a 
different isotopic combination of NH and/or Li. 

It is important to consider the overlap between the 
atomic and molecular traps. If both species are held in 



a magnetic quadrupolc trap and the atoms are initially 
much colder than the molecules, the atomic cloud will be 
much smaller than the molecular cloud. This may be ad- 
vantageous because it means that atom-molecule collisions 
will take place only at the low magnetic fields sampled by 
the atoms, and for many potentials such low-field colli- 
sions have low inelasticity even at temperatures around 
10 mK. However, there is a price to pay because poor 
overlap extends the time needed for sympathetic cooling. 
This might be alleviated by using a Ioffc-Pritchard trap 
to increase the size of the atomic cloud in the early stages 
of cooling. 



Acknowledgments 

We are grateful to Liesbeth Janssen for valuable discus- 
sions. This work is supported by EPSRC under the collab- 
orative project CoPoMol of the ESF EUROCORES Pro- 
gramme EuroQUAM. 



References 

1. M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wie- 
man, E.A. Cornell, Science 269(5221), 198 (1995) 

2. K.B. Davis, M.O. Mewes, M.R. Andrews, N.J. VanDruten, 
D.S. Durfee, D.M. Kurn, W. Ketterle, Phys. Rev. Lett. 
75(22), 3969 (1995) 

3. J. Doyle, B. Friedrich, R.V. Krems, F. Masnou-Seeuws, 
Eur. Phys. J. D 31(2), 149 (2004) 

4. R.V. Krems, Int. Rev. Phys. Chem. 24(1), 99 (2005) 

5. B. Friedrich, J.M. Doyle, Chem. Phys. Chem. 10(4), 604 
(2009) 

6. J.J. Hudson, B.E. Sauer, M.R. Tarbutt, E.A. Hinds, Phys. 
Rev. Lett. 89(2), 023003 (2002) 

7. H. Bethlem, W. Ubachs, Faraday Discussions 142, 25 
(2009) 

8. D. DeMille, Phys. Rev. Lett. 88(6), 067901 (2002) 

9. M.T. Bell, T.P. Softley, Mol. Phys. 107(2), 99 (2009) 

10. W.D. Phillips, Rev. Mod. Phys. 70(3), 721 (1998) 

11. H.L. Bethlem, G. Meijer, Int. Rev. Phys. Chem. 22(1), 73 
(2003) 

12. J.M. Hutson, P. Soldan, Int. Rev. Phys. Chem. 25(4), 497 
(2006) 

13. J.G. Danzl, E. Haller, M. Gustavsson, M.J. Mark, R. Hart, 
N. Bouloufa, O. Dulieu, H. Ritsch, H.C. Nagerl, Science 
321, 1062 (2008) 

14. F. Lang, K. Winkler, C. Strauss, R. Grimm, J. Hecker Den- 
schlag, Phys. Rev. Lett. 101, 133005 (2008) 

15. K.K. Ni, S. Ospelkaus, M.H.G. de Miranda, A. Pe'er, 
B. Neyenhuis, J.J. Zirbel, S. Kotochigova, P.S. Julienne, 

D. S. Jin, J. Ye, Science 322, 231 (2008) 

16. J.G. Danzl, M.J. Mark, E. Haller, M. Gustavsson, R. Hart, 
J. Aldegunde, J.M. Hutson, H.C. Nagerl, Nature Phys. 
(2010) 

17. J.D. Weinstein, R. deCarvalho, T. Guillet, B. Friedrich, 
J.M. Doyle, Nature 395(6698), 148 (1998) 

18. D. Egorov, W.C. Campbell, B. Friedrich, S.E. Maxwell, 

E. Tsikata, L.D. van Buuren, J.M. Doyle, Eur. Phys. J. D 
31(2), 307 (2004) 



10 



Alisdair O. G. Wallis et al.: The prospects of sympathetic cooling of NH molecules with Li atoms 



19. W.C. Campbell, E. Tsikata, H.I. Lu, L.D. van Buuren, 
J.M. Doyle, Phys. Rev. Lett. 98(21), 213001 (2007) 

20. E. Tsikata, W.C. Campbell, M.T. Hummon, H.I. Lu, J.M. 
Doyle, New J. Phys. 12(6), 065028 (2010) 

21. M.T. Hummon, W.C. Campbell, H.I. Lu, E. Tsikata, 
Y. Wang, J.M. Doyle, Phys. Rev. A 78(5), 050702 (2008) 

22. S. Hoekstra, M. Metsala, P.C. Zieger, L. Scharfenberg, J.J. 
Gilijamse, G. Meijer, S.Y.T. van de Meerakker, Phys. Rev. 
A 76(6), 063408 (2007) 

23. S.Y.T. van de Meerakker, R.T. Jongma, H.L. Bethlem, 
G. Meijer, Phys. Rev. A 64(4), 041401(R) (2001) 

24. D.J. Haxton, S.A. Wrathmall, H.J. Lewandowski, C.H. 
Greene, Phys. Rev. A 80, 022708 (2009) 

25. R. deCarvalho, J. Kim, J.D. Weinstein, J.M. Doyle, 
B. Friedrich, T. Guillet, D. Patterson, Eur. Phys. J. D 
7, 289 (1999) 

26. D.J. Larson, J.C. Bergquist, J.J. Bollinger, W.M. Itano, 
D.J. Wineland, Phys. Rev. Lett. 57(1), 70 (1986) 

27. A. Ostendorf, C.B. Zhang, M.A. Wilson, D. Offenberg, 
B. Roth, S. Schiller, Phys. Rev. Lett. 97(24), 243005 
(2006) 

28. C.J. Myatt, E.A. Burt, RW. Christ, E.A. Cornell, C.E. 
Wieman, Phys. Rev. Lett. 78(4), 586 (1997) 

29. A.G. Truscott, K.E. Strecker, W.I. McAlexander, G.B. 
Partridge, R.G. Hulet, Science 291(5513), 2570 (2001) 

30. G. Modugno, G. Ferrari, G. Roati, R.J. Brecha, A. Simoni, 
M. Inguscio, Science 294(5545), 1320 (2001) 

31. P. Soldan, J.M. Hutson, Phys. Rev. Lett. 92(16), 163202 
(2004) 

32. M. Lara, J.L. Bohn, D. Potter, P. Soldan, J.M. Hutson, 
Phys. Rev. Lett. 97, 183201 (2006) 

33. M. Lara, J.L. Bohn, D.E. Potter, P. Soldan, J.M. Hutson, 
Phys. Rev. A 75, 012704 (2007) 

34. M. Tacconi, E. Bodo, F.A. Gianturco, Theor. Chem. Acc. 
117, 649 (2007) 

35. M. Tacconi, E. Bodo, F.A. Gianturco, Phys. Rev. A 75, 
012708 (2007) 

36. M. Tacconi, L. Gonzalez-Sanchez, E. Bodo, F.A. Gi- 
anturco, Phys. Rev. A 76, 032702 (2007) 

37. P.S. Zuchowski, J.M. Hutson, Phys. Rev. A 78, 022701 
(2008) 

38. P.S. Zuchowski, J.M. Hutson, Phys. Rev. A 79, 062708 
(2009) 

39. P. Soldan, P.S. Zuchowski, J.M. Hutson, Faraday Discuss. 
142, 191 (2009) 

40. A.O.G. Wallis, J.M. Hutson, Phys. Rev. Lett. 103, 183201 
(2009) 

41. T.E. Mehlstaubler, K. Moldenhauer, M. Riedmann, N. Re- 
hbein, J. Friebe, E.M. Rasel, W. Ertmer, Phys. Rev. A 
77(2), 021402(R) (2008) 

42. P.J. Knowles, C. Hampel, H.J. Werner, J. Chem. Phys. 
99(7), 5219 (1993) 

43. J. Cizek, J. Chem. Phys. 45(11), 4256 (1966) 

44. T.H. Dunning, Jr., J. Chem. Phys. 90, 1007 (1989) 

45. H.J. Werner, P.J. Knowles, R. Lindh, M. Schiitz, P. Celani, 
T. Korona, F.R. Manby, G. Rauhut, R.D. Amos, A. Bern- 
hardsson et al., MOLPRO, version 2006.1: A package of 
ab initio programs (2006), see http://www.molpro.net 

46. S.F. Boys, F. Bernardi, Mol. Phys. 19(4), 553 (1970) 

47. C.R. Brazier, R.S. Ram, P.F. Bernath, J. Mol. Spectrosc. 
120(2), 381 (1986) 

48. T.S. Ho, H. Rabitz, J. Chem. Phys. 104(7), 2584 (1996) 



49. P. Soldan, J.M. Hutson, J. Chem. Phys. 112(9), 4415 
(2000) 

50. K.M. Jones, E. Tiesinga, P.D. Lett, P.S. Julienne, Rev. 
Mod. Phys. 78, 483 (2006) 

51. G.F. Gribakin, V.V. Flambaum, Phys. Rev. A 48, 546 
(1993) 

52. M. Mizushima, Theory of Rotating Diatomic Molecules 
(Wiley, New York, 1975) 

53. M.L. Gonzalez-Martinez, J.M. Hutson, Phys. Rev. A 75, 
022702 (2007) 

54. R.V. Krems, A. Dalgarno, J. Chem. Phys. 120(5), 2296 
(2004) 

55. H. Cybulski, R.V. Krems, H.R. Sadeghpour, A. Dalgarno, 
J. Klos, G.C. Groenenboom, A. van der Avoird, D. Zgid, 
G. Chalasihski, J. Chem. Phys. 122(9), 094307 (2005) 

56. W.C. Campbell, T.V. Tscherbul, H.L Lu, E. Tsikata, R.V. 
Krems, J.M. Doyle, Phys. Rev. Lett. 102(1), 013003 (2009) 

57. J.M. Hutson, S. Green, MOLSCAT computer program, 
version 14, distributed by Collaborative Computational 
Project No. 6 of the UK Engineering and Physical Sci- 
ences Research Council (1994) 

58. M.H. Alexander, D.E. Manolopoulos, J. Chem. Phys. 86, 
2044 (1987) 

59. L.M.C. Janssen, P.S. Zuchowski, A. van der Avoird, G.C. 
Groenenboom, J.M. Hutson, arXiv:1012.2814 (2010) 

60. L.M.C. Janssen, P.S. Zuchowski, A. van der Avoird, J.M. 
Hutson, G.C. Groenenboom, arXiv:1012.0804 (2010) 

61. P.S. Zuchowski, J.M. Hutson, Phys. Chem. Chem. Phys. 
13, 3669 (2011) 

62. E. Majorana, Nuovo Cimento 9, 43 (1932) 

63. W. Petrich, M.H. Anderson, J.R. Ensher, E.A. Cornell, 
Phys. Rev. Lett. 74, 3352 (1995) 



